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Abstract 

In conventional nuclear experiments a beam of accelerated nuclei collides 
with a target nucleus that is surrounded by other nuclei in a molecule, in 
condensed matter, or in a plasma environment. It is shown that for low 
collision energies possible nuclear reactions (including deuterium fusion) are 
strongly boosted by the environment. The effect originates from a chain of 
preliminary three elastic collisions which transform the projectile-target ex- 
periment into one with colliding beams. Firstly, the projectile-target pair of 
nuclei undergo elastic scattering in which the projectile shares its energy and 
momentum with the target nucleus. Then the projectile and target nuclei 
collide with different heavy nuclei from the environment. These latter colli- 
sions change the velocities of the target and projectile nuclei and set them 
again on the collision course. Finally, the same pair of nuclei collide inelas- 
tically, this time giving rise to the nuclear reaction. The increased energy of 
the target nucleus increases the relative velocity up to y/2 times, resulting in 
a drastic exponential increase of the probability to penetrate the Coulomb 
barrier and therefore sharply increasing the likelihood of the nuclear reaction. 
Applications to laser-induced fusion are discussed. 
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I. INTRODUCTION 

We suggest a new mechanism that strongly increases the probability of low-energy nuclear 
reactions in condensed matter, in molecules, or in a plasma environment. Although we do not 
aim to interpret a particular experiment, our interest in the subject was definitely inspired by 
a few publications, see |1|-§J and references therein, that claim a substantial increase of the 
DD fusion cross-section in solids. One obvious possibility for a boost of the fusion originates 
from the motion of the target deuterium caused by its vibrations in a condensed matter (or 
other) environment, as discussed in ||. Moreover, the vibration can be stimulated by the 
Coulomb interaction of the projectile deuterium with the target one, giving an additional 
increase for the probability of the fusion, as was found in ||. In this paper we suggest 
another effect which proves to be very effective in the low-energy region. 

Consider some nuclear reaction initiated by the collision of a pair of nuclei in a con- 
ventional beam-type nuclear experiment. Suppose that the target nucleus is surrounded by 
other heavy nuclei. This happens when it is deposited in a molecule, in condensed matter, 
or in a plasma environment. The energy of the projectile is supposed to be lower than the 
Coulomb barrier, and we are interested in the probability of a nuclear reaction. Remember 
that the Coulomb barrier makes the probability P of the reaction depend exponentially on 
the relative velocity v 12 of the colliding pair || 

where Z%, Z2 are the charges of the projectile and target nuclei, vn is their relative velocity. 
(Here and afterwards the indexes 1 and 2 mark the velocities, momenta, and energies of 
the projectile and target nuclei, respectively.) If we assume that the energy of the target 
nucleus E% is negligible, E2 ~ 0, then the collision velocity depends only on the energy of 
the projectile Ex, v 12 — V\ = J2Ei/m\, making the probability (|1.1|) depend exponentially 
on this energy. Let us imagine now that there exists some mechanism which can split the 
same amount of energy between the two colliding nuclei. Moreover, let us assume that this 
mechanism can also put the nuclei in a head on collision. Then obviously the relative velocity 
may be larger. Its possible maximum for a given energy E\ is v m = \jlE\j where /x is the 
reduced mass for the pair of nuclei. For the collision of two deuterons m\ = m 2 = 2/i, then 
the maximum possible relative velocity v m = \/2vi is substantially larger than the initial 
velocity, v m > v\. This increase of the relative velocity and a corresponding gain in the 
energy of the relative motion is due to reduction of the energy of the center of mass motion 
which was a half of the total energy in the initial projectile deuteron problem. Eq.( |l.l|) 
shows that collision with the velocity v m = \/2v\ makes the nuclear reaction much more 
probable. 

This simple observation inspires one to look for a mechanism which can fulfill two func- 
tions. Firstly, it should force the projectile to share its energy with the target. Secondly, it 
has to set the pair onto the collision course. The first goal can be achieved if we consider an 
elastic scattering of the projectile on the target nucleus. This collision obviously transfers 
part of the projectile energy to the target and can substantially enlarge the relative velocity 
of the two nuclei. This event should be considered as a preliminary collision, since we are 
seeking the possibility for the nuclear reaction to occur in the following inelastic collision. 
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(1.1) 



However, the necessary inelastic event cannot happen by itself. After the elastic collision, 
the two nuclei move in opposite directions in the center-of-mass reference frame, and there- 
fore cannot collide again, at least this is impossible in the vacuum. The situation changes 
when the environment is present. The important point is that scattering of the considered 
colliding particles on nuclei of the environment can change the directions of their propaga- 
tion. To simplify the consideration, let us assume that the nuclei of the environment are 
much heavier than the colliding pair (which is usually true for the most interesting case of 
light colliding nuclei). Then, the scattering on nuclei of the environment does not change 
the energies of the target and projectile nuclei, but changes the directions of their velocities. 
These new velocities may put the pair on the collision course. There arises therefore a possi- 
bility for the second collision of the same projectile-target pair. A sketch in Fig.|l] illustrates 
this idea. Since this final collision happens with larger relative velocity, it can initiate the 
nuclear reaction with much higher probability. 

The picture described needs a chain of three elastic collisions (TEC). The initial elastic 
collision of the pair accelerates the target, then the two collisions of the target and projectile 
nuclei with different nuclei of the environment change the relative velocity of the pair in such 
a way that the final inelastic collision becomes possible. One can find some resemblance of 
this process with the cannon (carambole) shot in the billiards game. 

The collisions of the projectile and target nuclei with atoms of the environment can, 
generally speaking, result in ionization of these atoms. One can expect, however, that an 
energy loss due to ionization is much lower than the collision energy, provided the latter is 
higher than the atomic ionization potential. That is why even if the ionization happens, 
the nuclei can still move along the trajectories which lead to their final inelastic collision. 
In other words, it is sufficient for our purposes to assume that collisions with atoms A and 
B happen quasielastically. The colliding nuclei do not lose a substantial amount of their 
energy, but transfer a large momentum. 

Proposing the TEC mechanism, we find that it greatly, by orders of magnitude, increases 
the probability of the nuclear reaction for a low projectile energy. In Section we find the 

we 



wave function for the colliding pair which takes into account TEC. In Sections |TJ and IV 



derive the amplitude and the probability of the nuclear reaction initiated by TEC. Section 
[V| examines the role of screening, while in Section [VI| we present a numerical example which 
shows how effective the TEC mechanism is for deuterium fusion. The maximum effect 
is achieved when the target deuterium is located at the center of the square formed by 
heavy atoms and the beam of projectile deuterons is parallel to the surface of this square. 
Conclusions are presented in Section |VII| . 



II. WAVE FUNCTION 

Our goal is to find the wave function which takes into account TEC. Note firstly that 
each step of TEC happens in some region, and the higher the collision energy, the more 
localized is the region. After each collision the active nucleus propagates a long distance 
before colliding again. Fig.|l] illustrates this statement. This indicates that each of the TEC 
collisions may be described in terms of the corresponding scattering amplitudes (physical 
amplitudes, i.e. on mass shell amplitudes). Let us call f p t{0 P t) the scattering amplitude 
for the elastic collision of the projectile with the target nucleus, and f p A(& P A), /ts(^ts) 
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the scattering amplitudes for collisions of the projectile nucleus with an atom A and the 
target nucleus with an atom B. We expect that the wave function of the target-projectile 
pair of nuclei which takes TEC into account S^ TECS) ip(r p , r t ) should be proportional to these 
amplitudes, 5 {TEC) ip(r p ,r t ) oc f pt {9 pt ) fpA(d pA ) /*b(0*b)- 

Furthermore, one can describe the propagation of the target and projectile nuclei in 
between the collisions in terms of the semiclassical approximation which is valid since the 
nuclei are heavy in comparison with electrons. Therefore we need to find the classical trajec- 
tories which describe this propagation. This many-body problem can be greatly simplified if 
we assume that in between the elastic collisions nuclei propagate along straight lines. This 
requires that the collision energy satisfies 

J?L^Z p -Z t ■ 27 eV, (2.1) 

llTLi 

where pi, m p are the momentum and mass of the projectile nucleus, and Z p , Z t are the charges 
of the projectile and target nuclei. Condition ( J2.1| ) is valid for any reasonable energy of the 
projectile. The relevant classical trajectories for both projectile and target nuclei start at 
some initial point r^, where the first projectile-target elastic collision takes place. From this 
point the projectile propagates with constant velocity v' x to the point where an atom A is 
located, and then, after scattering, acquires the velocity v'/ and moves to the final point r p . 
Similarly, the target moves from the point r« with the velocity \' 2 towards an atom B located 
at r# and then with the velocity Vg to the final point r t . Fig. ||] shows the Feynman diagram 
which illustrates this discussion. Obviously, the initial projectile-target collision preserves 
the total energy and momentum of the pair, while collisions of the projectile and target nuclei 
with atoms A and B preserve the absolute values of their momenta, provided the atoms A 
and B are sufficiently heavy. [] These classical conditions fix the classical trajectories for 
the projectile and target nuclei and allow one to express all relevant parameters, such as 
Yi , Vp, V t , Vp, v" etc, as functions of the final coordinates r p , r t and initial momenta pi, P2 
of the projectile and target nuclei. From the classical trajectories one can construct the 
corresponding semiclassical wave function 

^V P1 , P2 (r,r) = Fexp^, (2.2) 
F = X fpt(Qpt) fpA{0 P A) ftsiOtB) , (2.3) 

where S is the classical action for the described trajectories, and F is the corresponding 
pre-exponential factor. Eq. ( |2.3|) takes into account the fact that the wave function should 
be proportional to the scattering amplitudes for TEC. From the simple dimensional analysis 
we conclude that the only unknown parameter X in (|2.3|) should have the dimension of 
(length) -3 . The only available parameter of length in the problem is a separation of atoms 



Atoms A and B can be ionized during collisions. This process, generally speaking, reduces 
the total energy of the pair. However, this energy loss is comparable with the atomic ionization 
potential and, consequently, is much lower than the total energy of the pair. This allows us to 
neglect the ionization process in the following consideration. 
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a in the condensed matter (or molecule, or plasma) environment. Therefore we find an 
estimate X ~ 1/a 3 , which agrees with the more involved calculations discussed below. 

Remember that we need the wave function in order to describe the final inelastic collision 
of the pair. This requires that the projectile and target nuclei eventually come to the same 
point r, where the final collision takes place. Therefore it is sufficient for us to consider 
the wave function only at this point, i.e. when r p = r t = r. The method outlined above is 
discussed in detail for this specific case in Appendix A. The found wave function ( |A8| ) has 
the form 

<5 (T£C V Pa P2 (r, r) = MM ^4 exp % -S , (2.4) 

R(r) |r — ta\ |r — r^l n 



where the classical action S is given by ( |A9| ) . The wave function ( |2.4| ) is proportional to the 
scattering amplitudes, as was anticipated; the corresponding scattering angles 9 pt , 9 p a, QtB 
can be found from the classical trajectories. The wave function is also inversely proportional 
to the factors |r — r^|, |r — re|, as it should be since the waves that describe the projectile 
and target nuclei at the final point r originate from the points and r#, respectively. The 
factor R(r) in the denominator is an effective distance describing the propagation of the 
projectile and target nuclei from their common initial point, to and r#, respectively. 
Its explicit form is presented in Eq. ( |A10| ), but for our purposes it is sufficient to keep in 
mind that it is comparable with atomic separation 

R(r) ~ tab ■ (2.5) 

The semiclassical wave function ( |2.4| ) agrees with a simple estimate ( |2.2| ), ( |2.3| ). 

Eq. (|2.4j) is valid for the case when both the projectile and target nuclei in the initial 
state are described by plane waves. The validity of the plane wave approximation for the 
projectile nucleus is justified by its sufficiently high energy (|2.1| ). In contrast, the motion 
of the target nucleus arises due to vibrations of this nucleus in a molecule or in condensed 
matter, or by temperature reasons in a plasma environment. Therefore, generally speaking, 
the wave function of the target nucleus is not a plane wave. However, typical target momenta 
are small compared with the projectile momentum, 

P2 < Pi ■ (2.6) 

This fact allows one to modify the approach outlined above to accommodate the vibration 
motion of the target nucleus. Suppose that the target nucleus is described by the wave 
function 0t(r) in the initial state. Consider an event in which the target nucleus possesses 
a momentum p 2 before the collision. The amplitude of this event equals the the Fourier 
component of the target wave function, <p t (p 2 ). Multiplying ( |2.4| ) by this amplitude, we find 
the contribution which this event gives to the wave function of the pair. Integrating over all 
possible momenta of the target, we obtain the following wave function for the pair, 

^)* (ri ,r 2 ) = / ^V Pl , P2 (r,r)0> 2 )^ , (2.7) 

that takes into account a distribution of the target momentum in the initial state. This 
approach, justified by inequality (|2.6j ), is similar to the momentum approximation that is 
well known in nuclear physics @. 
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The semiclassical nature of the problem allows one to make the next step, fulfilling 



integration in (2/7) explicitly. Note that according to (|2.6| ) only small momenta p 2 con- 
tribute to the integral. Therefore we can consider the limit p 2 — > in the wave function 
S^ TEC ^ip pi )P2 (r, r). In this limit the classical action S , (r,p 1 ,p 2 ) should be expanded in the 
Taylor series 

S(t, pi, p 2 ) = S(r, pi, 0) + r, ■ p 2 + . . . , (2.8) 

where it is taken into account that the derivative of the classical action over the initial 
momentum gives the initial coordinate for the considered trajectory, 

8S(r a P " P2) = r, . (2.9) 
<9p 2 



Higher derivatives, denoted by dots in (|2.8|) , are negligible. The action is the only quantity 
in the wave function 5^ TEC ^il> Pl! p 2 (r, r) which needs to be considered accurately, using the 
series expansion. For other parameters one can simply take the limit p 2 = 0. As a result 
we derive from 



5™* PllP2 (r,r) =^V Pl ,p 2 =o(r,r)exp^(r l -p 2 ) . (2.10) 



Substituting this into (|2.7|) we find the final expression for the wave function 

^)*(r,r) = 5( TEC ^ pijP2= o(r,T)Un) , (2-11) 
where <pt{ Y i) is the wave function of the target nucleus calculated at the point r^, the latter 



one is a function of r as is discussed above (see also Appendix A). Eq. ( 2.11 ) has a clear 
physical meaning. In order to set the TEC mechanism in action, the initial collision of the 
projectile and target nuclei should happen at some particular initial point r,, starting from 
which the projectile and target propagate and undergo rescattering, arriving eventually at 
the final point r. Therefore, the wave function of the pair should be proportional to the 
amplitude to find the target located at the point r^. This later amplitude is simply equal to 
the wave function of the target nucleus (fit( r i)- Since the momentum of the target is much 



smaller than that of the projectile, we can neglect in (|2.11| ) the initial slow motion of the 



target nucleus g, and describe what happens after the initial collision with the help of the 
function 6^ TEC ^ip pitP2=0 (r, r). 

Eqs. ( p. 11| ), ( f2.4| ) are the main result of this section. They present the wave function 
for the TEC problem in terms of classical trajectories that describe the propagation of the 
projectile and target nuclei between the collisions, the scattering amplitudes responsible for 
the elastic collisions, and the wave function of the target nucleus. 



The role of this motion is investigated in some detail in ||. 
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III. AMPLITUDE OF NUCLEAR REACTION 



Let us call A(pi, p 2 ) the amplitude of the nuclear reaction initiated by a collision of the 
pair of nuclei with momenta pi and P2 in a vacuum. According to fll.lQ this amplitude 
strongly depends on the collision velocity. It is convenient to account for this dependence 
by introducing the amplitude £?(pi,p 2 ), 

( ite 2 ZiZ 2 \ 

4(Pi,Pa) = exp I — IB(pi,p 2 ). (3.1) 

Here the exponential factor describes the under-barrier penetration through the Coulomb 
barrier (that strongly depends on the velocity v = \pi/m p — p 2 /m t \), while £>(p 1 ,p 2 ) de- 
scribes the proper nuclear processes which happen afterwards. 

Let us express the amplitudes of events which happen in the environment (in condensed 
matter, molecules, etc) in terms of the vacuum amplitude. Consider first the direct collision 
of the projectile with the target, when one neglects preliminary elastic rescattering. In this 
case the wave function of the pair is simply a product of the incident plane wave describing 
the projectile nucleus and the target wave function, 

^ Pl (r p ,r t ) = exp -(p p ■ r p ) <p t (r t ) . (3.2) 

The amplitude of the direct nuclear reaction, in which the reaction products in the final 
state have momentum Q, can be written as 

M ( DIR \Q) = A(p u Q - pi) J ^ Pl (r, r) exp ~(Q ■ r) d 3 r (3.3) 
= A(pi,Q-pi)0 t (Q-pi) . 

Here the integration over the coordinate with the weight exp— -|(Q ■ r) ensures that the 
total momentum of the products of the nuclear reaction is equal to Q. 

Let us now take the effect of TEC into account. We need firstly to use the wave function 
( p.ll ) which describes the amplitude of the event in which the inelastic collision happens 



at the point r after TEC. Secondly, we need to multiply this amplitude by the amplitude 
of the nuclear reaction A(p", p 2 ') and, additionally, by the factor exp— |Q • r (that again 
ensures that the products of the nuclear reaction have the total momentum Q). Integrating 
the obtained result over the coordinate, we derive the following expression for the matrix 
element, 

5 {TEC) M(Q) = [ A(p''p")8^ TEC ^(r,r) exp--(Q.r) d 3 r . (3.4) 
J n 

The procedure used to derive this result is quite similar to the case of the direct collision 



( P-3| ). The only important distinction is the wave function, which for the case ( |3.4j ) accounts 
for TEC. The momenta p", p 2 of the projectile and target nuclei arise due to TEC and, 
henceforth, differ significantly from the initial momenta pi, p 2 . This fact drastically changes 
the exponential factor in ( p.l|) , as is discussed below. 
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IV. PROBABILITY OF NUCLEAR REACTION 



To calculate the total probability of the nuclear reaction, one needs to fulfill integration 
of the square of the amplitude over all possible momenta of the nuclear reaction products 
in the final state. For the direct collision this gives 



oc J \M^ Dm \Q) 



2 d 3 Q 



^(pi,Q-pi)<MQ-p 



2 d 3 Q 



(4.1) 



|A(pi,0) 



Deriving the last identity we took into account the fact that according to (|2.6|) the Fourier 
component of the target wave function is localized in a region of small momenta. Eq. ( |4.1| ) 
indicates that the probability of the nuclear reaction initiated by the direct collision in the 
environment is equal to the probability of the same event in the vacuum. 
Consider now the probability of the nuclear reaction initiated by TEC 



W/™oc J \6^M(Q)\ 



d 3 Q 



cxp 



27re 2 Z X Z 2 
hv'{ 2 



#(Pl,P2 



A(p'l,p^5^ E ^(v,v) d 3 r 

fpt(9pt) fpA^pA) /ts(^ts) , 



(4.2) 



R(r) 



r A 



r B 



in, 



d 3 r 



Here the first identity takes into account (|3.4|) rewriting the integration over the momenta 
d 3 Q as integration over the point of the final inelastic collision d 3 r. The second identity 
presents the wave function S^ TEC ^(r, r) explicitly using Q2.ll ), (|2.4j ), and (|3.1| ). The velocity 
in the exponential factor, 
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Pi 



n" 
Pi 

m t 



(4.3) 



depends on the momenta p", p 2 of the projectile and target nuclei after TEC. 



Let us show that the major contribution to the integral over d 3 r in 
the vicinity of a segment which connects and r B . Consider instead of 
relevant integral 



1 



exp 



2vre 2 Z X Z 2 



Uv'l 



d 3 r . 



originates from 
2f) an important 



(4.4) 
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Let us choose cylindrical coordinates with their origin at and the axis pointing along the 
vector Tba- If we ca U z the distance from the origin in the direction of this axis and p the 
radius vector in the orthogonal plane, then the radius vector can be written as r = (p,z). 
The projectile and target nuclei after TEC move along the segment r — and r — r^, 
respectively. Therefore the velocities in the cylindrical coordinate can be written as 



(4.5) 



(p, z tab) 
z 2 + {tab - z) 2 



S 



Expanding in powers of p/tab, one presents the collision velocity as a function of the coor- 
dinates, 



<2 = K + O 



1 

1 - - 



v[v' 2 



r 2 AB P 2 



2 (v[ + v> 2 ) 2 z* (r AB - zf 



(4.6) 



We assume in ( }4.5| ) that z belongs to the segment < z < tab- In this region the collision 
velocity is large, therefore its contribution is significant. Outside this region, when either 
z < or z > tab, the collision velocity is small, drastically reducing the integrand in ( |4.4| ) . 
Substituting (fOl) in ([4.4|), one finds 



rAB 



dz J cPpexp 



2vre 2 Z X Z 2 



hv t 



r A\ 



r B\ 



exp 



r 2 AB P 2 



z 2 {r A B ~ zf 



where v t is the sum of velocities, 



v t = v'l + v' 2 ' = v[ + v' 2 



(4.7) 



(4.8) 



The last identity here is valid because we assume that the collisions with atoms A and B 
are elastic ones (or quasielastic ones, if we have in mind that ionization of heavy atoms can 
occur). The dimensionless parameter £, defined as 



e 2 ZxZ 2 v[v' 2 



h 



3 ' 



(4.9) 



governs integration over d 2 p in ( f4.7|) . We consider low collision energies, therefore £ can be 
assumed to be so large that 



?r£ > 1 . 



(4.10) 



This inequality guarantees that separations from the z axis which are essential in the inte- 
grand in (|4.7|) are small compared to the separation between nuclei A and B, 



P^r A B ■ 



(4.11) 



All essential parameters which govern the wave function (including v[,v' 2 , etc) depend, gen- 
erally speaking, on the final point. However, in the vicinity of the z axis their dependence on 
p becomes insignificant, the only important dependence of the integrand on p is the factor 
p 2 in the exponent. This greatly simplifies integration over d 2 p in (|4.7|) , which results in 



1 



ir A B 



exp 



2yre 2 Z X Z 2 



dz . 



(4.12) 



A similar method applied to (^2) gives 

// tfsi2 \ fpt{Qpt) fpA{9pA) fpB{QpB) 



W^)oc/L4(p»p* 



B?{v)t\ b 



dz 



(4.13) 



Here the momenta are defined for the head-on collision which happens on the segment 
connecting and r B , 
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Pi = "Pi— ' P2 = P2— • (4-14) 

tab tab 

The scattering amplitudes in (|4.13|) are also defined for the head-on collision. 

|2 



Notice that the probability ( |4.13| ) is proportional to the density |0t(rj)| of the target 
nucleus at the point of initial collision. To make the probability bigger, we need, therefore, 
to distinguish the case when the target nucleus is well localized in the condensed matter 
or molecular environment. In this, most interesting for us, situation, further simplification 
of ( |4 . 1 3| j becomes possible. When the target nucleus is localized we should suppose that 



the relevant classical trajectories, which describe propagation of the projectile and target 
nuclei, start from some point in the vicinity of the localization point. Our previous results 
indicate that the trajectories finish at some point on the segment tab- These facts altogether 
fix the trajectories for the projectile and target nuclei uniquely. As a result, the angles in 
the scattering amplitudes as well as all velocities v' l ,v' 2 ,v t in the integrand in (|4.13|) can 
be considered as parameters defined for the given classical trajectories. In this case, the 
probability of the nuclear reaction initiated by TEC W TEC , Eq. ( }4.13| ), can be presented as 



W^oc/CKpMI 2 , (4-15) 

where the coefficient is 

1 | fpt {Opt) fpA{6pA) fpB {QpB)\ 2 f i , / \ 1 2 



K= f^>~~ /W*. (4.16) 



The factor /C, that appears in (|4.15|) and is defined in (|4.16| ), is the most important result 
of the calculations above. It can be conveniently rewritten in terms of the differential cross 
sections which describe TEC, 

We introduce here the density of the target nucleus n t (ri) = |0i(r.;)| 2 at the point of the 
initial elastic collision rv In order to fulfill integration in (|4.17| ) over the z coordinate of 



the point of the final collision, one needs to remember that rj is a function of the point of 
the final collision r, see discussion in Section [II] and Appendix A. This integration cannot 
be simplified any further in the general case, but it presents no difficulties for numerical 



estimations, see example in Section [VI[ The consideration above ignored the temperature of 
the target. If we assume that the temperature is lower than the projectile energy, T <C E\, 
then by slightly modifying the approach above, we arrive at the same result ( |4.17| ), with the 
density n t (r-j) a function of the temperature. 

Eq. ( }4.17| ) shows that each scattering process gives a factor into the probability which 
is identical to the elastic cross section; that result could have been anticipated without 
any calculations. Having established this, one could believe that an estimate for the pre- 
exponential factor should look like /C ~ ^m~(@pt) ~^t(^pa) ^f~(^tB)^/ r AB > because the 
atomic separation tab seems to be the only relevant parameter in the problem. However, 
our calculations (^4.17|) reveal that there exists an additional enhancing factor k, 

^wM^^^f (418) 

k = r 2 AB J nt(Ti)dz ~ r 2 AB n 2 t /3 (0) ~ ~ 100 , (4.19) 
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where nt(0) = 1/a 3 is the maximum density of the target nucleus at the point of localization, 
and a is the corresponding length which describes the region of localization. The latter one 
depends on the amplitude of vibrations exhibited by the target in the considered condensed 
matter or molecular environment. This length can usually be estimated as ~ 10% of an 
atomic separation ||, which justifies the last identity in ( }4.19| ). f] Eqs. ( |4.18| ), ( |4.19[ ) predict 
that 

K oc oc n 4 / 3 , (4.20) 

where n is the density of the condensed matter environment. We will use this fact for 
estimations described below in Section |V11| . 

We consider above a chain of scattering processes with the participation of atoms A 
and B, indicating this fact in the notation Kab in ( |4.17p . If the process happens in large 



molecules, or in condensed matter, there may be several heavy atoms located close to the 
target nuclei. Each of them provides the possibility for the necessary elastic rescattering. 
Therefore, if we want to consider all possible rescattering processes we need simply to add 
the contributions arising from scattering on different neighbor nuclei. If the nucleus A 
participates in the process then it gains a large momentum. Therefore, the contributions 
from scattering on different nuclei do not interfere, and we need to sum the probabilities. 
The result again can be presented with the help of ( |4.15| ), where the coefficient is 



IC = J2^ab. (4.21) 



A,B 



The summation here runs over all neighbor nuclei of the environment. 

Consider the ratio of the probability of the nuclear reaction initiated by TEC to the 
probability of the direct event, 

w( TEC ^ 

n = worn • ( 4 - 22 ) 

Calculating this ratio with the help of fl4.1p , ( (4.15[ ), one should remember that the proba- 
bility of the nuclear reaction includes the sophisticated matrix elements which describe the 
proper nuclear processes and are defined as the amplitude -B(pi, P2), Eq. (|3~1"|). These latter 
factors vary on the scale of nuclear energies which are much higher than the collision energy. 
Therefore, with high accuracy these nuclear matrix elements can be considered as energy- 
independent constants, which cancel out in the ratio ( |4.22[ ). Having this fact in mind, we 
find the following expression for the ratio 1Z, 



1Z = JC exp 



2vre 2 ZxZ 2 ( 1 1 



(4.23) 



Here the factor JC is defined in ( |4.17| ), ( |4.21| ). Note an important inequality v m > v%, which 
guarantees that the ratio 1Z increases exponentially for low energies E\. 



Note that in the derivation of ( |4.17 ), the length of localization a was assumed to be sufficiently 

dU ' dQ. 



large, a 2 S> da £^ , ^ff-, which is true for sufficiently high energy. 
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The applicability of the developed theory is restricted by several bounds. One of them 
comes from within the theory. Inequality ( 4.10 ) requires that £ is large, thus putting a bound 
for the energy from above. For an important example of the collision of two deuterons, this 
inequality gives V\ < (e 2 /h) 7r/(4\/2), resulting in E\ < 15 keV. From below, the energy in 
our approach is obviously restricted by a typical potential energy in a crystal or molecular 
environment, that necessitates (12.11). 



V. SCREENING 

Screening plays an essential role in low energy collisions of nuclei in condensed matter 
@|. Let us summaries briefly relevant facts. When the two colliding nuclei come sufficiently 
close, their influence on the electrons of the environment can be described by their total 
Coulomb charge Z = Z^ + Zb- This charge for low energy collision creates discrete atomic- 
like energy levels for electrons. The electrons of the environment can populate these levels 
and form an effective atom around the combined nucleus. This is the strongest possible 
response which the electrons can exhibit in screening the Coulomb field. An occupation 
of discrete levels in the combined atom results in energy production. The maximal energy 
produced in this process is equal to 

U = E B (A Z )-Z- A , (5.1) 

where Eb{Az) is the total binding energy of the atom Az that has Z electrons surrounding 
the combined nucleus with nuclear charge Z, and A is the binding energy for each electron 
in the given condensed matter environment. The energy variation U can be transferred 
into the kinetic energy of the colliding nuclei. This transfer may be expected to be efficient 
for small collision velocity. Eq. ( |5.1|) shows that screening makes the collision energy higher 
E co i — > E co i + U. This, in turn, makes the relative velocity larger 

via -> v s = ^v 2 l2 + 2U/n , (5.2) 

thus increasing the probability to penetrate through the Coulomb barrier ( |1.1|) and boosting 
the nuclear event. Screening in this description is described by the screening potential 
U. The larger U, the stronger is the screening and the higher is the probability of the 
nuclear reaction. It is clear that the given consideration may overestimate a role of screening 
because, firstly, a probability to populate the ground state of the effective atom should, 
generally speaking, be less than unity, and, secondly, because the energy production should 
not necessarily be transferred solemnly into the kinetic energy of the colliding nuclei, but 
may be used for excitation of the environment as well. These two factors could make 
screening potential less than what is predicted by (|5.1| ). Accurate account of all relevant 
factors is difficult and, probably, due to this reason Ref. || uses U as a parameter to fit the 
experimental data. 

The above discussion shows that screening results in an increase of the relative velocity. 
This effect takes place in the direct collision as well as in the TEC process. To take screening 
into account we need to increase velocities in both these collisions. Applying this procedure 
to the relative probability ( |4.23| ) of the nuclear process initiated by TEC we find 
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1Z — > IZs = /C exp 



2vre 2 ZiZs 
ft 




(5.3) 



where vi t s and ^5 are defined in accord with ( |5.2j ) 



Vl,S = 



\lv\ + 2U/fi , 



(5.4) 
(5.5) 



It is clear that screening produces stronger impact for the direct collision because the relative 
velocity in this case is smaller than in the collision initialed by TEC. Therefore screening 
should reduce the relative probability IZs (|5.3j ) in comparison with its non-screened value 
1Z ( |4.23j ), see Section [VT| for more details. 

VI. APPLICATIONS FOR DEUTERIUM-DEUTERIUM FUSION 

Let us apply the results obtained above to the important case of deuterium fusion d+d — > 
3 H + p in a collision of the projectile deuterons with target deuterons implanted in a solid. 
Eq. ( [4.17| ) shows that we need to estimate the cross sections for TEC. For the initial deuteron- 
deuteron elastic collision one can use the conventional Rutherford formulae 



where E\ is the energy of the projectile deuteron and 9 p t is the scattering angle in the 
center-of-mass reference frame. It should be noted that the target deuterium in a condensed 
matter environment can form a neutral atom; the originally charged projectile nuclei prop- 
agating through the environment can also pick up an electron becoming neutral as well. 
Thus, during the deuteron collision the charges of the nuclei are partly screened by the 
deuterons' electrons. This fact, however, does not change estimate ( |6.1|) which assumes the 
pure Coulomb interaction between the nuclei. The necessary significant variation of veloc- 
ities of the deuterons requires that their collision should be sufficiently hard, i.e. happen 
with small impact parameter b, b ~ e 2 /E p . Therefore, for energies higher than the atomic 
energy, E p 3> 27.21 eV, the impact parameter must be much smaller than the Bohr radius, 
b ^ aQ. The smallness of the impact parameter shows that screening produces a small effect 
for deuterium scattering. 

The same argument can be applied to the collision of a deuteron with a heavy atom. 
Electrons located at the deuteron nucleus or in an outer shell of the heavy atom cannot 
produce any significant effect because major events during the scattering happen for small 
separation between the deuteron nucleus and the nucleus of the heavy atom. Q 



4 This argument is in line with the discussion of the role of ionization during collisions, see Section 
||. The collision with the heavy atom can, generally speaking, be accompanied by ionization of 
the heavy atom or the deuteron atom. However, the ionization process involves an outer shell of 




(6.1) 
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The above arguments are valid for an outer shell of the heavy atom. In contrast, elec- 
trons of inner shells are localized closer to the nucleus and therefore produce a strong effect, 
providing screening for the Coulomb field created by the heavy nuclei for those small sep- 
arations which are responsible for the scattering. The screening can be conveniently taken 
into account in the semiclassical approximation, the validity of which is justified by the large 
mass and energy of a deuteron. In this approach one needs firstly to calculate the deuteron 
scattering angle as a function of the impact parameter p 







71 




2U{r) 



-1/2 



dr 



(6.2) 



Then the cross section is found from the conventional relation QIC 

— (f)) 9 dp 
dQ { ' sinOde ' 



(6.3) 



We perform these calculations using the well-known Thomas-Fermi atomic potential || for 
the atomic potential energy U(r) = Utf(j"). Fig. [3] presents results for deuteron scattering 
from gold atoms. Our calculations for other atoms give similar results. The cross section is 
comparable with an area with the atomic dimension Oq for low energies and monotonically 
decreases with energy. 

To proceed further we need to fix the geometry of the collision. Let us assume that the 
momentum of the projectile deuteron p x is orthogonal to the segment AB, i.e. pi • tab — 0. 
Assume also that the target deuteron is located in the vicinity of the point most convenient 
for this case, r = r^ B /2 — (p 1 /p 1 )r^ B /2. Then the scattering angles are 



e. 



P t 



90 c 



e 



pA 



IB 



135° 



(6.4) 



In this geometry, shown in Fig. f|, the final inelastic collision happens at the point r^#/2. To 
be specific, we assume that the atoms A and B are the atoms of gold (though this assumption 
is not essential, other heavy atoms would produce a similar effect). 

The velocities after the first deuteron-deuteron collision are v[ — v' 2 — Vi/y/2, which 
makes the collision velocity for the final inelastic collision 



From this we deduce for the ratio (|5 



1Z = /C exp 



v t 



2vre 2 



V2 




(6.5) 



(6.6) 



the heavy atom or the deuteron atom. This means that the ionization happens when the deuteron 
nucleus is separated by a large distance ~ ao from the heavy atom nucleus. In contrast, large 
variation of the velocity of the deuteron nucleus happens at small separations, where ionization 
does not play a role. We conclude again that scattering with ionization can be considered as a 
quasi-elastic process which does not affect significantly the collision with large variation of the 
deuteron velocity. 
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where 



/C — C — — 7- _n 

£r AB El 



d(JpA 

dVL 



(135°) 



(6.7) 



and C = 0.7 is a numerical factor. []. For the considered case one finds that the parameter 
£ defined by ( |4.9|) is equal to 

1 e 2 

£ = -^-^- . (6.8) 



The velocities u^g and u tjl g in the exponent of ( |6.6| ) are defined in (|5.4|) , (|5.5|) , thus taking 
screening into account. 

Let us estimate tab as a typical atomic separation, say the separation between atoms of 



gold metal, tab = 407.82 nm = 7.71 au. The enhancing coefficient k according to Q4.19 ) is 
k ~ 100. 

Fig. § shows the results of our calculations for the ratio 7Z ( |6.6|) of the probability of the 
nuclear reaction initiated by TEC to the probability of the direct collision. The coefficient 
/C was calculated from (|6.7Q as explained above. The screening potential U, which influences 
the velocities v^s and v t> s m C I5.4|) , ( ]5.5f) was found experimentally in Ref. || for the case 
of Au to be U = 22.8 ± 11.0 eV. Estimation (|5.2|) predicts a higher value. Note that the 
combined atom in the case considered is He, its binding energy is E(He(ls 2 )) = 79.0 eV. 
If we assume that A ~ 13.6 eV (binding energy of the Hydrogen atom), than ( |5.2| ) gives 
U = 51.8 eV. Deviation of this result from the experimental data prompts us to present in 
Fig. |5] data which show variation of the enhancing factor 1Z both on the collision energy and 
the screening potential U. We see that for collision energies below ~ 500 eV the ratio 1Z 
becomes large, 1Z > 1, showing that the TEC makes the probability of the nuclear reaction 
bigger than in the direct inelastic collision. For lower energies this enhancement becomes 
stronger. Screening makes the nuclear reaction more probable, but it results in the decrease 
of the relative probability probability 1Z for processes initiated by TEC. This damping effect 
arises because screening manifests itself stronger for the direct process, in which the collision 
energy is lower. The stronger screening, i.e. the higher the potential U is, the less important 
role plays the TEC mechanism. 

The analysis above assumes that the solid state environment keeps the target deuteron 
localized in the center of a square which is cornered by heavy atoms, and the beam of the 
projectile deuterons moves parallel to some side of this square. This geometry is the most 
favorable because it provides an enhancement factor (|4.18|) . For an arbitrary geometry, the 
enhancement ( 4.18| ) does not take place, making the probability of TEC lower by a factor 
of ~ 100. 

We conclude that the fusion below 500 eV is strongly influenced by the TEC mechanism. 



5 It is calculated as a product C = 2 • (2/3) • c in which the coefficient 2 takes into account the 
fact that the sum in ( 4.21| ) includes two possible trajectories for the pair of nuclei which lead each 
one of them towards either atom A or B, the coefficient 2/3 originates from Eq. ( |A12 ), and the 
numerical factor c = (y/2 + 1/2)" 1 = 0.522 arises from the integration over dz = cdzi in Eq. ( p7|) 
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VII. CONCLUSION 



We have demonstrated that there is a boost in the probability of a nuclear reaction, in 
a condensed matter or molecular environment, when three elastic collisions (TEC) precede 
the inelastic collision of two nuclei. When operative, the TEC mechanism provides an expo- 
nential enhancement for nuclear reactions with low collision-energy. One should remember, 
however, that this boost does not fully compensate for the suppression due to the Coulomb 
repulsion. It accounts, roughly speaking, for only 30% (0.3 ~ 1 — l/y/2, see (|6.6| ) ) of the 
exponent which describes tunneling through the Coulomb barrier, leaving 70% of this expo- 
nent intact. Moreover, the compensation originates from a chain of several (three) collisions 
which produces an additional pre-exponential, but nevertheless important, damping factor. 
Additionally, screening also reduces the role played by TEC. As a result of this we find that 
the probability for a nuclear reaction is enhanced by TEC only for sufficiently low collision 
energies. An example of DD fusion discussed in Section [VI] demonstrates that the effect of 
TEC manifests itself strongly for energies below 500 eV, being negligible for higher energies. 
Obviously, in this low energy region the probability of fusion is very low, by a factor of 10 14 
lower than for the energy 2 keV. This should make direct measurement of the discussed 
effect very difficult in the near future. [] 

Nevertheless, the effect results in a strong boost for low-energy nuclear reactions which 
may be important for different applications. One of them can be looked for in laser-induced 
fusion, in which ions can be accelerated by the laser field, see e.g. |§. This acceleration 
can, generally speaking, create a flux of deuterons which move from a surface inside a target 
exposed to a laser field. If these ions have energies below 0.5 keV, than the TEC mechanism 
can boost the fusion. One can deduce a practical conclusion from this. To allow the TEC 
mechanism to manifest itself in laser-induced fusion, the deuterium target should contain a 
sufficient amount of heavy atoms. 

One can look for other applications of the TEC mechanism, for example, in astrophysics. 



APPENDIX A: WAVE FUNCTION IN THE SEMICLASSICAL 

APPROXIMATION 



Let us derive the wave function which describes the projectile and target nuclei with 
TEC taken into account. This goal can be achieved with the help of the Feynman diagram 
presented in Fig. |2|. The solid lines with open arrows on the diagram show the projectile 
and target nuclei. Their propagation in the intermediate states should be described with the 
help of the corresponding Green functions G p , G t - Since we assume that the energies of the 
nuclei are sufficiently high, see fl2.1|, we can use the free-particle Green functions G a = G^\ 
a = p,t. In the coordinate-time picture the propagator of a free particle is 

^^(s^f^H 3 ^)- (A1) 



6 Hence, the discussed effect can not be directly responsible for some unexplained anomalies in 
fusion which have been observed experimentally. 
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Open circles on the diagram describe interactions which are responsible for TEC. Let us 
use notation M a b (with appropriate labels ab = pt, pA, tB) to describe the corresponding 
matrix elements. The important point is that each step of TEC happens inside some localized 
region, each collision being well separated from other ones. As a result, each step of TEC 
can be described by a physical ("on mass shell") scattering amplitude. The matrix elements 
M a b are related to the scattering amplitudes f a b by conventional formulae, 

f ab = -^M ab , (A2) 

where m ab is a reduced mass of the pair of particles ab. The scattering angles in these 
amplitudes are specified below. 

The considered diagram Fig. |2| gives the following contribution to the wave function, 

rtf rtf roo r 7 

5(™)*( ri) r 2 ) = / dt A dt B dt f / d 3 r i exp-(P-r i + M / ) (A3) 

J JO J max (t A ,t B ) J Tl 

x M pt Gf (r i} 0; v A , t A )M pA G p °\r A , t A ; r^^Gf (r,, 0; r B , t B ) M tB Gf (r B , tB', I"2, tf) 

rtf ft 



rtf rtf roo r ^ 

/ dt A / dt B / dtf / d 3 riF exp —S . 

JO Jo J max Ha ,tn) J n 



I max (tAja) 

The integration here includes integration over the point of initial collision and integration 
over the time variables. The latter describe the times of the elastic collisions: the projectile 
collision with the atom A (t a ), the time of the target collision with atom B (t B ), and the 
final time (tf). The quantities P = p x + p 2 and E = E\ + E 2 are the total momentum of 
the pair in the initial state and the corresponding energy. 

We use in ( |A3| ) the "time-coordinate" presentation for the propagators which is conve- 
nient for our purposes since locations of atoms A and B are fixed. Another, probably more 
widely used, presentation provides the "energy-momentum" picture, in which propagations 
of particles are described by conventional energy denominators. We verified that (|A3|) can 
be directly derived from the energy-momentum presentation. 

The quantity S in the exponent in (|A3|) originates from the phases of the four Green 
functions (AT) combined with the phase factor P • r, + Etf, 



S = 5(r 1; r 2 ; t A , t B ,t f , v i )=P-r l + Etf (A4) 
m p \r A - r^| 2 m p \r - r A \ 2 m t \r B - r^ 2 m t \r - r B \ 2 

i(t A -uf + 2(t f -t A y + 2{t B -uy + 2(t f -t B y ' 

One recognizes in this expression the classical action which describes the free motion of the 
two nuclei between the collisions. The factor F in ( |A3| ) describes all relevant pre-exponential 
factors. 

The heavy, high energy nuclei have a wavelength much less than the atomic separation. 
Therefore, propagation of the nuclei in a condensed matter environment can be described 
semiclassically. This means that the integral in ( |A3| ) can be evaluated using the saddle-point 
approximation. The saddle-point position should be found from the conditions 

dS_ _ dS_ _ dS_ _ dS_ _ 

d^ i ~dT A ~dt B '~drf~ 1 [ ' 
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which coincide with the classical equations of motion for the two nuclei. The classical 
equations specify the locations and moments of time of TEC which are the functions of the 
coordinates ri, r 2 . They specify also the momenta of the projectile and target nuclei in the 
intermediate state. These momenta allow one to find the scattering angles that govern the 
matrix elements M a b in ( |A3| ) . 

In the semiclassical approximation, the final pre-exponential factor is related to the 
action, since 



J Fexp—S dt A dt B dtfd 3 ri 



2ttz\ 3/2 F 
IT/ v/del^ 77 



(A6) 



Here 5"' is the 6x6 matrix of second derivatives of the classical action over the six integration 
variables t A ,t B ,tf,Ti. To simplify the formulae, we restrict the following consideration to 
the specific case when ri = r 2 = r. In this case one finds that the matrix 5"' is 



S" 



(A + £ 





V 



2E[ 

tfA 
Ppa 
t A 





2E> 2 U + J- 

z \*S tfB 
_2Ei L 

tfB 
Pice 
t B 



2E[ 

tfA 

2E> 2 



2 



tfA 

( JEL _|_ E 'l 

\tfA tf B 




t A 

tB 


b ^ \JI + Tb-) J 



(A7) 



where a, (3 = 1,2,3 describe the Cartesian projections of the vectors. Calculating the 
determinant of the matrix ( |A7| ) and substituting the result in and ( |A3| ), we find the 
wave function which can be written in the following form 



6 TEC ^(r,r) 



fptjdpt) fpA^pA) fwiPtB) i Q 

R{r) |r — |r — r B \ h 



(A8) 



The action on the classical trajectory can be presented in a convenient, more familiar form 



S = P • r 

The factor -R(r-j) in 

R(v) = U 



p' x {\y - r A \ + \r A - Yi\) +p 2 (|r - r B | + |r B - r<|) . 

) comprises the result of the calculation for detS"', 

\ 1/2 



/ 2 

; 12 



+ 



m p t Bi + m t t Ai 



sin 2 7 



m p t Bi + m t t A j 
m„ + m t 



(A9) 



(A10) 



where the angle 7 is defined by 



cos 7 



vi-v 2 

v[v' 2 



(All) 



The notation above assumes the usual convention t ab = t a — tb- Eq. (|A10|) looks complicated, 
but it greatly simplifies for a particular important case, m p = m t , t A i = t B i, and cos 7 = 0, 
when it becomes 



R(r) 



-r AB . 



(A12) 
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Otherwise, in the general case, one can make an estimate R(r) ~ tab, which suffices for our 
purposes. 

The wave function fl^5|), that takes into account TEC, is the final result of this section. 
To avoid confusion, let us note that it does not account for the Coulomb repulsion of the 
nuclei when they finally collide inelastically. This effect is taken care of separately, see the 
corresponding suppression factor exp[— 27ie 2 ZiZ 2 / (hv" 2 )] in Eq. (|4.2|) . 
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FIGURES 



FIG. 1. The sketch illustrates the main idea of three elastic collisions. Lines with arrows 
show classical trajectories of the projectile and target nuclei that are depicted by small circles. 
The projectile nucleus collides firstly with the target nucleus, then the projectile and target nuclei 
are scattered by atoms A and B depicted by large circles. After that the inelastic collision of the 
projectile and target nuclei takes place, resulting in the nuclear reaction. The relative velocity in 
this latter collision is larger than the initial velocity of the projectile nucleus, which enhances the 
probability of the event. 

FIG. 2. The Feynman diagram which describes the mechanism of three elastic collisions. Solid 
lines with arrows in the intermediate states show propagators of the projectile and target nuclei. 
They are marked by the corresponding momenta. Open circles describe the three elastic scattering 
processes. Each scattering happens at some well defined point: the initial one takes place at r-j. 
The scattering of the nuclei by atoms A and B occurs at the points and r# where these atoms 
are located. The final inelastic collision takes place at rj = r. Since all these points are well 
separated in space, these collisions may be described by the physical scattering amplitudes (on 
mass-shell amplitudes). Calculation of this Feynman diagram is presented in Appendix A. 

FIG. 3. The differential cross section for elastic scattering of deuterons by gold atoms 

calculated with the help of Eq. ( |6.3| ) for the scattering angle 8 = 135°. 

FIG. 4. Three elastic collisions for deuteron fusion. The initial momentum of the projectile 
deuteron pi is orthogonal to the segment Rab- The target deuteron is located so that the final 
inelastic collision happens with velocity v\2 = V2v\. This is the largest possible collision velocity, 
therefore the considered geometry is the most suitable one. 

FIG. 5. The relative probability of the DD fusion initiated by TEC. The ratio 1Z, defined in 
(|4.22j ), is presented versus the projectile energy. Target deuterons are implanted in an environment 
containing gold atoms which are responsible for rescattering, the geometry of collisions coincides 
with the one shown in Fig. |j. Calculations are based on Eq. (|6.6| ) with the elastic cross section 
from Fig. ||. The screening potential: thick solid line - no screening U = 0, thin solid line - U = 10 
eV, dashed line - U = 20 eV, dashed-dotted line - U = 30 eV, dotted line - the strongest screening 
U = 50 eV. 
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